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ABSTRACT 

The dynamical friction problem is a long-standing dilemma about globular clusters 
(hereafter, GCs) belonging to dwarf galaxies. GCs are strongly affected by dynamical 
friction in dwarf galaxies, and are presumed to fall into the galactic center. But, GCs 
do exist in dwarf galaxies generally. A solution of the problem has been proposed. If 
dwarf galaxies have a core dark matter halo which has constant density distribution 
in its center, the effect of dynamical friction will be weakened considerably, and GCs 
should be able to survive beyond the age of the universe. Then, the solution argued 
that, in a cored dark halo, interaction between the halo and the GC constructs a new 
equilibrium state, in which a part of the halo rotates along with the GC (co-rotating 
state) . The equilibrium state can suppress the dynamical friction in the core region. In 
this study, I tested whether the solution is reasonable and reconsidered why a constant 
density, core halo suppresses dynamical friction, by means of N-body simulations. As 
a result, I conclude that the true mechanism of suppressed dynamical friction is not 
the co-rotating state, although a core halo can actually suppress dynamical friction 
on GCs significantly. 

Key words: methods: N-body simulations - galaxies: dwarf - galaxies: kinematics 
and dynamics - galaxies: star clusters - galaxies: structure. 



1 INTRODUCTION 

The world we live in is a hierarchical universe, in which 
galaxies are made by a myriad of merging events. A large- 
scale numerical simulation based on the Cold Dark Matter 
(CDM) theory has been operated as a greatly declarative 
method for the hierarchical scenario. It indicated that the 
structure of the universe develops from the small dark mat- 
ter clumps which collapsed first, and result in the formation 
of large and massive dark mater halos. In such a formation 
history, it is appropriate to consider that dwarf galaxies are 
fundamental 'building-blocks' and expected to be the old- 
est structures of the universe. Dwarfs are believed to have 
important clues in understanding the hierarchical universe. 

In this paper, I will discuss the dynamical friction 
problem which refers to orbital motion of GCs in dwarf 
galaxies. The drag force of dynamical friction is negligibly 
weak for GCs in the Milky Way. In contrast, it operates 
strongly in small system s like dwarf galaxies (see chap. 8 of 
iBinnv fc TremamelBoOsh . Thus, the GCs in dwarfs are pre- 
sumed to lose their orbital energy and fall into the galactic 
center by strong friction force from the dark matter halo. Ac- 
cording to results of both analytical and numerical studies, 



E-mail:inoue@astr. tohoku.ac.jp 



the timescale f or a GC to fall into the cen ter is of the or- 
der of ~ lGyr (Trcmaine ll976l : Hernandez 
Oh et aLlbOOOtrV cspcrim 2000l 



Gilmordll998l; 
^, .20011 : ISanchez-Salcedo et al.l 
20061 : iGoerdt et al.ll2006h . Nevertheless, even in the present 



universe, these GCs still do exist and keep their orbital mo- 
tions. For example, the Fornax dSph galaxy has five GCs 
which are metal poor and as ol d as the universe, thus r esem- 
bling the GC s of Milky Way (iBuonanno et al.l [l998l. Il999l; 
Strader et all 120031 : iMackev fc Gilmord liool iGreco et all 



2007) 



However, by an analytical approach, 

Hernande z fc Gilmord (| 1998T ) have discovered that a 
King model halo can significantly weaken the effect of 
dynamical friction in the core region. As for a cuspy halo 
(NFW profile or singular isothermal sphere), a GC is sucked 
into the galactic cente r by th e dynamical friction (see the 



fig. 2 of Goerdt et al. (|2006l) ) . The analytical approach 



of iHernandez fc GilmorehfT99a ) was constructed on the 



Chan drasekhar dynamical friction formula l|Chandrasekhar1 
Il943f). On t he o ther hand, by N-bo dy simulations, 
IGoerdt et all (I2006D and lRead et"atl (|2006l ) (hereafter, R06) 
confirmed the cessation of dynamical friction on a GC in 
a core region of halos. However, R06 concluded that an 
important key to this suppressed dynamical friction is a 
'co-rotating state'. They argued that a part of halo particles 
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in the constant density core begin to rotate with the GC. 
The authors suggested that this dynamical state is a new 
equilibrium state including the GC. The dynamical friction 
ceases under this equilibrium; he nce, the GC could su rvive 
beyond the age of the universe. iGoerdt et all (2006) and 
R06 confirmed that these results don't depend on the mass 
of a GC, the orbital parameters (circular or elliptical orbit) 
of a GC, or the size of core structure of a halo (core radius). 

However, the concl usions of these studies, 
iHernandez fc Gilmord (|l998l ) and R0 6, imply a discrepancy 
betwe en them. The approach of IHernandez fc Gilmorej 
(1998) is based on the Chandrasekhar formula; hence, 
it cannot take account of velocity anisotropy of the field 
particles, because the formula is based on the assumption 
of isotropic velocity state. But, contrary, R06 concluded 
that the mechanism of suppressed dynamical friction is the 
very anisotropy in the velocity state: the co-rotating state. 
My aim in this paper is to assess the authenticity of the 
co-rotating state proposed by R06 as the mechanism of the 
suppressed dynamical friction. 

Although the result of R06 appears to be convincing, 
the new equilibrium, the co-rotating state, will be vulnerable 
to perturbation on the GC. Once the GC orbit is perturbed 
and the orbital plane inclined, the co-rotating equilibrium 
state will be broken easily. This means that the system is no 
longer equilibrium: the dynamical friction force on the GC 
will be rejuvenated. R06 have also mentioned this fragile na- 
ture of the co-rotating equilibrium state. In R06 they studied 
single GC cases only. But, real dwarf galaxie s don't neces- 
sarily have only one GC , but several or mo re |Durrell et al.l 
1 1999 : iMiller et all 1 19981 ; lLotz et all 120041 ). The GCs in a 
dwarf will be perturbed by the other GCs. Then, I infer 
that some GCs would fall into the galactic center by such 
orbital perturbation, and may merge and form a stellar 
nucleus cluster at the galact i c cen ter (|Miocchi et al.l [20061 : 
ICapuzzo-Dolcetta fc Miocchil l200si ). Actually, some dwarf 
galaxies have a nucleus stellar cluster at the center, and 
some observational researchers have disc ussed that some of 
these nuclei ma y be remnants of GCs (|Miller et alj Il998l ; 
lLotz et aHl2004r ). I conducted N-body simulations to exam- 
ine whether the GCs in dwarfs fall into the galactic center by 
perturbation from the other GCs, even in the case of cored 
halo structure. 

This paper is organized as follows: in the next section, 
I will explain my simulation method and models in detail; 
in the 3rd section I will show my simulation results and 
their analysis; I will discuss my results, comparing them with 
preceding studies, and give my conclusion in the 4th section. 



2 THE SIMULATIONS 

My simulational settings are almost the same as the N-body 
simulation of R06. The simulations are pure N-body simula- 
tions (no gas component). I use a Barne s-Hut modified tree- 
code ((Barnes fc Hutlll986l ; lBarnej|l99d ) in order to lighten 
the heavy burden of gravitational force calculation, setting 
an open angle of 9 = 0.5. A special-purpose calculator for 
collisionless N-body simulations, GRAPE-7 model 600, is 
used with the tree algorithm to accelerate gravity calcula- 
tion (|Makinolll99lf ). The total number of timesteps is 11841 
for the whole of a simulated period which corresponds to 10 



Gyr in real timescale. It takes roughly a half month to finish 
each simulation. 



2.1 The setting of halo model 

To imitate R06, I adopt the same spherical density distribu- 
tion: 
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with a = 1.5, (5 = 3.0, 7 = 0.0. The density in the core po 
is O.WMqpc~ 3 . The scale radius r s is set to 0.91kpc. The 
density is nearly constant at the center within 200-300 pc, 
which defines the core region. The virial mass of the halo is 
2.0 x 10 9 Mq. I add an exponentially decaying envelope to 
prevent instability at the outer region caus ed by an artificial 
cut-off radius r„; r l|Springel fc Whitelll999h . Velocity disper- 
sion of particles is given by the solution of Jeans equations 
as a function of radius, 
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where f3 is the anisotropy parameter. Although the effect of 
dynamical friction is sensitive to the velocity distribution of 
the field particles, isotopic velocity state is supposed to be 
reaso nable in inner region of dwarf halos l|Mashchenko et all 
2008). In this paper, I assume the isotropic velocity state in 
the halo, setting f3 = (oy = ag — o^). With this assump- 
tion, Eq[2] reduces to 
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The velocity distribution is determined by the local 
Maxwellian approximation, 



F(v,r) = 47r 



IwJ " eXP l-2^J' 



(4) 



where F (y, r) is a prob ability distribution function of 
velocity l|Hernquistl Il993f ). EqH is normalized so that 

JT F (v, r)dv = 1 

Lik e R06 andlGoerdt et all (|2006f ). I adopt a three-shell 
model l|Zemp et al.l |200"8T ), which consists of finer grained 
particles in inner regions and coarser particles in outer re- 
gions. This technique enables it to reduce computational 
run-time, and resolve much smaller scales in the inner region. 
But, this multi-shell model inevitably admits the heavier 
particles coming from the outer shells into the inner shells, 
and may induce two-body relaxation between these differ- 
ent mass particles. To avoid such an unfavorable artificial 
effect, I refine the particles in the outer shells depending on 
their orbit. From an set of initial position and velocity, I can 
calculate the pericenter distance of a specific particle in the 
smooth potential given by the density profile, Eq[T] By the 
pericenter distance, I detect the heavier particles which are 
supposed to intrude into the inner region. I divide these in- 
truding particles into a set of particles which have the same 
mass resolution as the particles in the inner region. The new 
particles will have the same radial velocity component as the 
original particle but a new random tangential component of 
the same magnitude as the original one. The divided par- 
ticles are randomly placed on a sphere whose radius is the 
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same as the initial galactocentric distance of the original par- 
ticle. For a simple explanation, let's suppose that the halo 
consists of shells A, B and C, from inner to outer shell. The 
shells A, B and C consists of particles which have the mass 
of rriA, tub and mc, respectively (mj < mj < mc). For 
example, if a certain particle which has the mass of mc and 
an initial position in the shell C intrudes into the shell A, the 
particle will be divided into a set of new particles which have 
the mass of tua, and the number of the ne w particles will 
be mc/tnA (for detail, see lZemp et alJlboOfif ): therefore, not 
all particles in a outer shell have a uniform mass resolution, 
but two-body relaxation can be minim ized in the inner mos t 
region. In the simulations of R06 or iGoerdt et al.l (2006), 
they have been missing the particle-dividing step; therefore, 
at this point my simulations have an improvement over the 
preceding studies. 

Specifically, the particle masses are raj = I7.8M0, 
raj = 356M0 and mc = 7118Mq. The inner most region, 
the shell A, is within 300pc (to be accurate, the shell A is not 
a shell, but a sphere). The shell B, the middle shell, is the 
region from 300 — llOOpc. The shell C, the outer most shell, 
is the entire region of outside of the shell B. The number of 
particles is 8.31 x 10 6 in total. The softening lengths of the 
particles, mj, and mc, are 3, 8 and 22 pc, respectively. 
I checked that my results were not sensitive to these values. 

The three-shell model has coarse resolution in the outer 
shells. Such heavyer particles may affect the nature of dy- 
namical friction. For confirmation, I ran a simulation with 
an uniform mass particle model, and checked that the results 
were not sensitive to these mo del settings (see Appendix A) . 

By analytical calculation, ISanchez-Salcedo et al.l (|2006T ) 
has proposed that the dynamical friction induced by the 
stellar component is not negligible in a cored dwarf halo. 
But, the purpose of this paper is to judge the co-rotating 
state proposed by R06 as the mechanism of the suppressed 
dynamical friction. For fair comparison with the simulation 
of R06 in which stellar component was excluded, I don't 
take the effect from stellar component into account in my 
simulations here. 



2.2 The setting of globular clusters 

In my simulations, each GC is represented by a point 
mass with mcc = 2.0 x 10 M®. The softening length 
is lOpc. Just to make sure, I ran another simulation in 
which a GC was resolved by many particles, and con- 
firmed that tidal disruption didn't destroy the GC. In this 
study, I don't consider mass-loss from a GC, merging be- 
tween GCs, or dynamical heating by halo potential. Al- 
though these effects may play important roles in the case 
of the resolved GCs jFuiii et al.l 120061: iMiocchi et al.ll2006l ; 
lEsauivel fc Fuchsll2007l ; ICapuzzo-Dolcetta fc Miocchfoosl ) . 
I consider the G Cs as point-masses f or the sake of compari- 
son with R06 or lGoerdt et~ail i|2006u . 
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Figure 1. The comparison between analytical and numerical re- 
sults of time evolution of orbital radius of a GC. The smooth 
lines indicate the analytical results, the waving jaggy lines in- 
dicate the N-body results. I estimated the analytical dynamical 
friction with In A = 3.0. The initial orbit of the GC is a circular 
orbit from radius 600 pc or 1 kpc. 
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3 THE RESULTS 

In this paper, I operate simulations for the cases of 1, 5 and 
30 GCs. In single GC cases, I give the GC a set of specific 
orbital parameters in each case. On the other hand, in multi 



Figure 2. The results of single GC cases, each line indicates 
different simulation. The red line represents the case in which a 
GC is placed on the circular orbit of radius 600 pc initially. The 
green, blue lines are for the GC initially placed on elliptical and 
radial orbit from radius 600 pc, respectively. The core region is 
inside 200-300 pc, in all simulations here. 



GC cases, I set their orbits at basically random as detailed 
in the following subsections. 



3.1 The single GC cases 

To begin with, I conduct some single GC cases. Initial orbit 
of the GC is circular and set to 600 pc or 1 kpc. In Fig[T] I 
show two comparison cases with the Chandrasekhar dynam- 
ical friction formula (the derivation of the analytical result 
is described in Appendix B .). As seen in the figure, the 
analytical results correspond to the N-body results before 
the GC enters into the core region (inside 200-300 pc). But, 
when entering into the core, these results diverge abruptly: 
the analytic results continue to fall into the galactic center, 
while the orbital shrinkage by dynamical friction stops in 
the N-body cases. 
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Figure 3. The density profiles in the case of the circular orbit 
with 1 kpc initial orbital radius (the upper jaggy lino in Fig[Tjl. 
The solid line is the initial state, and the dashed line indicates 
the end of the simulation (after 10 Gyr). The horizontal axis is 
distance from the galactic center, the ordinate is mass density of 
the halo, p[M Q /pc 3 ]. 

However, as I noted above, the effect of dynamical fric- 
tion is sensitive to the velocity distribution of the field par- 
ticles. But I confirmed that the core stalling of dynamical 
friction was not sensitive to details of the velocity distribu- 
tion function with various anisotropy parameters /3 in Eq(2] 

Next, I investigate the relation of the suppressed dy- 
namical friction with orbital eccentricities. The initial orbit 
is circular, elliptical or radial and set to 600 pc. In the case of 
elliptical orbit, the rotational velocity of the GC is initially 
set to 0.6v c (v c is the circular velocity at the initial position) 
and the radial velocity is 0. In the case of radial orbit, the 
GC is at rest initially. The results are indicated in Fig[2] 
As shown in the figure, when the GC enters into the core 
region, the orbital shrinkage by dynamical friction stops in 
all cases regardless of their initial orbital eccentricities. Af- 
ter the cessation of dynamical friction, the orbit expands a 
little again. This phenomenon is called 'kickback effect', and 
a detailed investiga tion about the effect has been done by 
iGoerdt et ail (120081 ). In this study, I don't treat this effect. 
These behaviors of the GC in the cored profile (the cessation 
of orbita l shrinkage, the kick back effect) are consistent with 
R06 and lGoerdt et all i|2006T ). The mass included in the halo 
core is heavier than a GC by two ord ers, which should be 
enough to operate dynamical friction. IGoerdt et alj (|2006l ) 
has confirmed independence of suppressed dynamical fric- 
tion from the size of the core region, r„. Moreover, I find 
that the density profile of the halo is scarcely changed by 
the GC (Fig|3} • The energy conservation rate of the system, 
1-E end /E ini , is 6.46 x TO" 3 . 

This result confirms the suppressed dynamical friction 
on a GC in a cored halo. However, what is the cause of 
it? In R06, the authors have proposed that a part of halo 
particles in core region are made to rotate with the GC by 
gravitational interaction (see the fig. 4 in R06). They called 
this dynamical state 'co-rotating state' and concluded it to 
be the mechanism to suppress the dynamical friction on the 
GC. I will examine whether the co-rotating state is the true 
mechanism or not. To identify the co-rotating state, I follow 
the same manner as R06. In order to visualize the velocity 



state of the halo, I sieve out the particles which have the 
radial distance r p — 150 — 300pc , subject to a condition, 



where J P:C means the specific angular momentum of a halo 
particle and the GC, respectively. A criterion parameter 6 
is set to 10°, which is the same value as R06. Because the 
halo potential is spherically symmetric and the GC is signif- 
icantly heavier than any halo particles, the direction of J c is 
assumed to be constant in time. The condition, Eq(5] screens 
out the field particles for which the direction of angular mo- 
mentum vector coincides with that of the GC within 8. Fig|J] 
indicates histograms of rotational velocity distributions in 
the case of a circular orbit for which the initial orbital ra- 
dius is 1 kpc (the upper jaggy line in Fig[T|. The upper 
panel of Fig(4] shows the initial state, the bottom panel is 
for t = 8.2Gyr (after the cessation of dynamical friction). As 
seen in the bottom panel, the co-rotating state is constructed 
after the dynamical friction is suppressed. The velocity dis- 
tribution becomes somewhat anisotropic: the fraction of pro- 
grade rotating particles seems to increase, whereas retro- 
grade particles decrease. The figure is consistent with the re- 
sult of R06. The over-plotted dashed lines in the histograms 
represent Gaussian fitting given by the minimum \ 2 method. 
In the bottom panel (co-rotating state), the peak height of 
the fitting line for pro-grade side h+ and the retrograde side 
h- is 0.0132 and 0.0106, respectively. The residual fraction 
of these, {h+ — h-)/h-, is 0.245. This value means that the 
peak of the pro-grade side is 24.5 per cent higher than that 
of the retrograde side. The minimum value of \x\ for the 
bottom panel is 0.00175, and |x|/ft,_ is 0.165. 

From this analysis, the existence of co-rotating state 
seems to be confirmed. But it may be marginal because the 
value of [h+ — h-)/h- is not significantly larger than \x\/h~. 
One point to note is that the direction of J c , which is as- 
sumed to be constant, actually fluctuates due to the N-body 
nature of this simulation (i.e., a finite number of particles 
has been used) For the sake of more precise discussion, I 
evaluate (h+ — h-)/h- and for slightly different di- 

rections of J c . I re-analyze the velocity states, changing the 
inclination of the vector J c in Eq[5] little by little. By this 
procedure, I draw contour maps of the value of (h+~h- ) /h- 
and \x\/h-. Fig[S] indicates the results. As shown in the fig- 
ure, the co-rotating state ((h+ — h-)/h- > 0) can be found 
within ~ 10° from the maximal direction. But, the value of 
(h+ — h-)/h- should be compared with the value of \x\/h- . 
The contour map of is given in the bottom panel of 

Fig(5] From the map, it is found that the range of the value 
is 0.12 < < 0.2 for the entire area plotted. There- 

fore, the value of (h+ — h-)/h- is almost submerged under 
the fitting error, \x\/h- > i n most of the area except the cen- 
tral region. This discussion means that the co-rotating state 
shown in the upper panel of Fig[5] is statistically marginal 
and unimportant, maybe except the central region (~ 3° 
from the maximal direction). Such weak anisotropy couldn't 
be expected to affect the dynamical friction on a GC. 

Finally, I inspect the co-rotating state for the depen- 
dence on radial distance from the galactic center. So far, 
I've analyzed the co-rotating state in the radial range of 
150pc < r < 300pc. I additionally carry out the same analy- 
sis in other radial ranges, r < 150pc, 300pc < r < 450pc and 
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Figure 4. The histograms of the distribution of azimuthal ve- 
locity for the halo particles satisfying the condition of Eq[5] 
and 150pc < r p < 300pc. The vertical axis is mass fraction. The 
upper panel is for the initial state, and the bottom panel is for 
t = 8.2Gyr. The positive side of represents pro-grade rotating 
particles with the GC, the negative side represents retro-grade 
motion. The dashed lines are Gaussian fitting for each side. 



450pc < r < 600pc, with J c unchanged. Fig[5] indicates the 
results. From the figure, it can be seen that the co-rotating 
velocity state is constructed only in 150pc < r < 300pc. 



3.2 The 5 GCs case 

From the discussion of the previous subsection, the co- 
rotating state not be influential in dynamical friction. Even 
if the co-rotating state suppresses the dynamical friction on 
a GC, because it would be fragile against orbital perturba- 
tion on the GC as proposed in R06, the dynamical friction 
would be rejuvenated by the presence of the perturbation. 
Actu ally, because real d warf galaxies generally have some 
GCs l|Miller et al.lll998l ). there are probably frequent per- 
turbations on the GCs in real dwarfs. As an additional test 
for the authenticity of the co-rotating state, I perform the 
simulations of multi GC cases. 

In this subsection, I present the result of the 5GCs case. 
The GCs are represented as point masses. The initial po- 
sitions of the GCs are randomly determined in the radial 
range of 300pc < r < 2kpc. The energy of each GC is also 
given randomly in the range of E < <3>2fcp C with a random 
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Figure 5. The contour maps of (/i_|_ — h-)/h- (upper) and 
|x|//i— (bottom). The center is the direction for which the value 
of (/i_|_ — h-)/h— is the largest ('the maximal direction'). The 
dotted circles indicate 5°, 10° and 20° from this direction, from 
inside to outside. 
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Figure 6. The radial dependence of (H+ — h—)/h—. The width 
of each bin represents a radial range of the analysis. J c is fixed 
to the maximal direction. 



direction of the velocity vector ($2fe pc is the potential energy 
of the halo at r = 2.0fcpc). To facilitate direct comparison 
with single GC cases, one GC (reference GC) is placed on a 
circular orbit with radius 600 pc. 

The results for all GCs are shown in Fig[7] The top 
left panel indicates the reference GC, and the other panels 
represent the other 4 GCs which have randomly chosen posi- 
tions and velocities. As shown in Fig[7] all orbits of the GCs 
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Figure 7. The simulation result of the 5 GCs case, showing the variations of orbital radius. The top left panel 
represents the reference GC which is initially placed on a circular orbit with a radius 600 pc. The others indicate 
GCs whose initial distance and velocity are randomly given. 



are frequently perturbed. From the comparison with single 
GC cases (FigHJ), the perturbation is expected to be caused 
by mutual interaction among the GCs, because the differ- 
ence between these simulations is the number of GCs only. 
Despite these perturbations, no GCs fall into the galactic 
center, with all GCs surviving and keeping their orbital mo- 
tions. This means that the dynamical friction is suppressed 
in a core region even under frequent perturbation. This im- 
plies the inconsistency with the description in R06: the vul- 
nerability of the co-rotatig state to perturbations. This re- 
sult casts a doubt on the co-rotating state as the mechanism 
of suppressed dynamical friction, together with the analysis 
of the single GC case. Actually, the co- rotating state can not 
be found in velocity histograms of the 5GCs case. 



3.3 The 30 GCs case 

To confirm the result of the 5GC case, I investigate the case 
of 30GCs. The number of GCs in this case is somewhat 
too large, because the act ual number of GC s in any dwarf 
galaxies is ~ 20 at most (|Miller et alj 1998). The settings 
for the initial condition is basically the same as the 5GCs 
case: random positions and velocities are given except for 
one GC which is placed on a circular orbit with radius 600 
pc (reference GC). Some simulation results of these GCs are 
shown in Fig [8] 

As seen from Fig(8] these results are essentially the same 



as 5GCs case. Dynamical friction is suppressed in the core re- 
gion despite more frequent perturbation than the 5GC case. 
Moreover, the co-rotation state is not confirmed in this case 
also, and can be guessed to be broken by perturbation among 
numerous GCs. The result reinforces my argument. 



4 DISCUSSION & CONCLUSION 

My simulation results indicate as follows. On one hand, dy- 
namical friction is indeed suppressed in a constant density 
core region, orbital shrinkage of a GC stops and the or- 
bital motion is sustained. But, on the other hand, this re- 
sult doesn't depend on the number of GCs. This means that 
the dynamical friction doesn't work in a core, even if the 
co-rotating state is broken by frequent perturbation. More- 
over, in single GC cases the co-rotating state seems to be too 
marginal to affect the dynamical friction. My main conclu- 
sion here is that the co- rotating state is not the true mecha- 
nism. There may be another reason why dynamical friction 
ceases in a constant density halo. 

However, I confirmed that a cored halo structure cer- 
tainly can weaken the dynamical friction force on a GC in 
the core region. This means that the cored structure can 
be the solution of the dynamical friction problem. Actu- 
ally, as an observational fact, the rotation curves of Low 
Surface Brightness (LSB) galaxies seem to be like a solid 
body rotation which means a nearly constant density distri- 
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Figure 8. The same as Fig[7] but for the 30GCs case. The top left panel represents the reference GC. 



bution, although CDM cosmological N-body simulations in- 
dicat e cuspy density distr ibutions l|Dalcanton fe Bernstein! 
l200d ; Ide Bl'ok et al] l200lf). Moreov e r, by the cosmological 
SPH simulation. iMashchenko et al.l (I2008I ) has discovered a 
dwarf galaxy which has a constant density core in the halo. 
They concluded that massive stars inject large amount of en- 
ergy into the dark matter halo via supernova explosions, and 
the feedback induced by bursty star formation can turn the 
cuspy density distribution of the dark matter into a cored 
profile in the halo center. The suppressed dynamical friction 
would reinforce the existence of a core structure in dwarfs. 
If it's the case that the actual dwarfs have a core region in 
the center, all GCs belonging to a dwarf are included within 
the core region of the dwarf. Thus, it can be expected that 
in a dwarf the largest galactocentric distance of the GCs in- 
dicates the minimum value of the core radius. For instance, 
in the Fornax dSph, the furthest GC fr om the galactic cen- 
ter h as a projected distance of ~ l.&kpc |Mackev fc Gilmorel 
2003). The core region of the Fornax would be larger than 
~ l.Gkpc. 

The most fundamental examination of dynamical fric- 
tion is e mbodied by the Chan drasekhar dynamical friction 
formula (|Chandrasekhad 1 1943ft . However, several approxi- 
mations have been dared in the derivation. Current astro- 
physicists realized that the formula is not always correct 
because of complex nonlinear effects(lJiang &: Bmnevll2000l ; 
lHashimoto et al.l I2003I ; iFuiii et al.l I2006I ; iKim et all l200Sft . 
Probing the dark matter distribution of dwarfs requires more 
sophisticated approaches to the nature of dynamical friction. 



ACKNOWLEDGMENTS 

The numerical simulations reported here were carried out 
on GRAPE systems kindly made available by CfCA (Cen- 
ter for Computational Astrophysics) at the National As- 



tronomical Observatory of Japan. The numerical code was 
based on a software distributed on the website of Joshua 
E. Barnes (http://ifa.hawaii.edu/~barnes/software.html). I 
thank Masafumi Noguchi for his fascinating and helpful dis- 
cussion. 



REFERENCES 

Barnes J. 1990, J. Comput. Phys., 87, 161 

Barnes J., Hut P. 1986, Nature, 324, 446 

Binny J., Tremaine S., 2008, Galactic Dynamics Second 

Edition. Princeton Univ. Press, Princeton 
Buonanno R., Corsi C. E., Zinn R., Pecci F. F., Hardy E., 

Suntzeff N. B., 1998, ApJ, 501, 33 
Buonanno R., Corsi C. E., Castellani M., Marconi O, Pecci 

F. F., Zinn R., 1999, AJ, 118, 1671 
Chandrasekhar S., 1943, ApJ, 97,255 
Capuzzo-Dolcetta R., Miocchi P., 2008, ApJ, 681, 1136 
Dalcanton J. J., Bernstein R. A., 2000, AJ, 120, 203 
de Blok W. J. G., McGaugh S. S., Bosma A., Rubin V. C, 

2001, ApJ, 552, 23 
Durrell P. R., Harris W. E., Geisler D., Pudritz R. E., 1996, 

AJ, 112, 972 
Esquivel O, Fuchs B., 2007, MNRAS, 378, 1191 
Fujii M., Funato Y., Makino J., 2006, PASJ, 616, 642 
Goerdt T., Moore B., Read J. I., Stadel J., Zemp M., 2006, 

MNRAS, 368, 1073 
Goerdt T., Read J. I., Moore B., Stadel J., 2008, preprint 

(astro-ph/0806.1951) 
Greco C. et al., 2007, ApJ, 670, 332 

Hashimoto Y., Funato Y., Makino J., 2003, ApJ, 582, 196 
Hernandez X., Gilmore G., 1998, MNRAS, 297, 517 
Hernquist L., 1993, ApJ, 86, 389 
Jiang I.-G., Binney J., 2000, MNRAS, 314, 468 



8 S. Inoue 




Figure Al. The comparison of the uniform mass model with the 
three-shell model. The result of the three shell model is the same 
as the case of the circular orbit from 1 kpc in Sec. 3.1. 



Kim H., Kim W.-T., Sanchez-Salcedo F. J., 2008, ApJ, 
679,33 

Lotz J. M., Miller B. W., Ferguson H. C, 2004, ApJ, 613, 
262 

Mackey A. D., Gilmore G. F., 2003, MNRAS, 340,175 
Makino J. 1991, PASJ, 43, 621 

Mashchenko S., Wadsley J., Couchman H. M. P., 2008, Sci, 
319, 174 

Miller B. W., Lotz J. M., Ferguson H. C, Stiavelli M., 

Whitmore B. C, 1998, ApJ, 508, 133 
Miocchi P., Capuzzo-Dolcetta R., Di Matteo P., Vicari A., 

2006, ApJ, 644, 940 
Oh K. S., Lin D. N. C, Richer H. B., 2000, ApJ, 531, 727 
Read J. I., Geordt T., Moore B., Pontzen A. P., Stadel J., 

2006, MNRAS, 373, 1451 
Sanchez-Salcedo F. J., Reyes-Iturbide J., Hernandez X., 

2006, MNRAS, 370, 1829 
Springel V., White D. M. 1999, MNRAS, 307,162 
Strader J., Brodie J. P., Forbes D. A., Beasley M. A., 

Huchra J. P., 2003, AJ, 125, 1291 
Tremaine S., 1976, ApJ, 203, 345 
Vesperini E., 2000, MNRAS, 318, 841 
Vesperini E., 2001, MNRAS, 322, 247 
Zemp M., Moore B., Stadel J., Carollo M., Madau P., 2008, 

MNRAS, 386, 1543 



orbit and the stalling of dynamical friction in the core re- 
gion. 



APPENDIX B: THE DERIVATION OF THE 
ANALYTICAL DYNAMICAL FRICTION 

The Chandrasekhar dynamical friction formula is 

dv c {r) _ 4tt In AG 2 p{r) M GC \ r ,„ s 2X 



dt 



erf(X) 



(Bl) 



where X = v c /(V2a(r)). Mgc is the GC mass. v c (r) is 
circular velocity of the GC as a function of radius. G is the 
gravitational constant. A = G(jm+biQc) ' wnere ^™ ls the 
largest impact parameter, m is a mass of field particles, Vo is 
a typical veloci ty of the system. In t his pa per, I approximate 
lnA = 3.0 fsee lBinnv fe Tremaine! |2008t) ). 

I approximate that the orbit is circular. Therefore, L = 
Mgc tv c and 



dL dr 
_ =Mg cv c -. 



(B2) 



Moreover, because dynamical friction drag force is parallel 
to v c , and therefore a torque 



dL , , v c 

IE = MGcr df 

From these equations, 

dr r v c 
dt v c dt 



(B3) 



(B4) 



Substituting Eci lBll for Eci lB4l I calculate the time evolution 
of a GC orbit by numerical integration of Eq |B4l 



APPENDIX A: THE REASONABILITY OF THE 
THREE-SHELL MODEL 

The three-shell model I used here is a complex model. To 
make sure that the heavier particles with large force sofening 
lengths in the outer shells don't affect the dynamical fric- 
tion unfavorably, I run a simulation with an uniform mass 
resolution model. The uniform mass model consists of par- 
ticle rriB described in Sec. 2.1. The total number of particles 
is 6.57 x 10 6 . The velocity state of the model is isotropic 
(Eel}. 

The result is shown in Fig lAll As seen, there is no dif- 
ference between the two models in time evolution of the GC 



